Multitime controlled linear PDE 

systems 



Cristian Ghiu and Constantin Udriste 

University Politehnica of Bucharest, Faculty of Applied Sciences, 
Department of Mathematics II, Splaiul Independent ei 313, 060042 
Bucharest, Romania, e-mail: crisghiu@yahoo.com 

University Politehnica of Bucharest, Faculty of Applied Sciences, 
Department of Mathematics-Informatics I, Splaiul Independent^ 313, 
060042 Bucharest, Romania, e-mail: udriste@mathem.pub.ro, 
anet . udri@yahoo. com 

Abstract 

We derive new results regarding the controllability and the reacha- 
bility of multitime controlled linear PDE systems of first order. These 
systems describe some important multitime evolution in engineering, 
economics and biology. Some of them come from evolution PDEs 
of superior order. The original results include a refinement and a 
supplement of multitime optimal control theory, developed in some 
recent papers by the second author. They refer to the complete in- 
tegrability conditions, conditions for the existence of solutions, path 
independent curvilinear integrals, the multitime fundamental matrix, 
multitime adjoint Cauchy problems, control space, controllability and 
reachability of phases, controllability gramian, reachability gramian, 
controllability matrix, counter-examples and commentaries. 
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1 Introduction 

Here a controlled system is a dynamic multitime linear PDE system on which 
one can act by using appropriate controls. Among the most common prob- 
lems that appear when studying such systems are multitime controllability 
problem and multitime reachability problem. 
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The multitime controllability refers to pairs of states that can be moved 
from the first one to the second one and the multitime reachability operates 
on the reverse order of states. Of course, the order of states is given by 
the product order (partial order) on multitime source space. The study of 
controllability of dynamical systems represented by normal PDEs starts in 
the papers [J2] - [23] , [2] (multitime maximum principle) , [7] , [SJ , [24] (max- 
imum principle in the context of weak derivatives), [4] (numerical methods 
for robust control), [TO] . [TT] (reachability of hybrid systems). Though many 
of situations are rather well understood, there are still quite challenging open 
problems due to the fact that the product order relation on multitime space 
is not total. 

This paper deals with control theory for systems governed by multitime 
linear PDE systems (m- flows). Section 2 presents a new and complete frame- 
work for the multitime nonautonomous linear PDE systems of first order. 
Section 3 contains original results about controllability and reachability of 
the controlled multitime nonautonomous linear PDE systems of first order. 
The controllability and the reachability of multitime autonomous linear PDE 
systems of first order is analyzed in Section 4. The comments (Section 5) 
show that in some other situations can occur strange mathematical phenom- 
ena due to the discontinuity of controls in multitime evolutions. 



2 Nonautonomous linear PDE system of 
first order 

We start with some mathematical ingredients related to evolution PDEs (m- 
flows). Let t = (t 1 ,. ••>*"*) G R m , called multitime, x = (x 1 ,...,x n ) T G 
R n = M ni i(R), and G C R m x ]R n be an open subset. We consider the 
evolution PDE system 



— (t)=X a (t,x(t)), Va = l,m, (2.1) 
where X a : G -> R n = M n ,i(R), X a = (X' a , . . . , X™). 

Definition 2.1. The PDE system (12. ip is called completely integrable if 
V(t ,x ) G G, 3D C R m , D open with t G D and 3x : D -> R n , x dif- 
ferentiate, such that (t,x(t)) G G, Wt G D , x verifies (12.11) and x(t ) = x . 

The following Theorems, 12.11 to I2.4[ represent new versions of some well- 
known results 0, [T21 [25]. 
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Theorem 2.1. Suppose the components X a are of class C l , Va = l,m. 

i) Any solution of the PDE system fl2.ll) is of class C 2 . 

ii) If the PDE system (12.11) is completely integrable, then 

^(t, x) + z)§§^, x) = ^(t, x) + X'(t, x), (2 2) 

V(t, a;) G G, Va, — l,m 

or in matrix notations 

dx a f dx a dx a \ _ dx^ (dx^ dxA 



d& ydx 1 dx n J p dt a v ' ' \dx x dx n J ' (2.3) 
V(t, x) G G, Va, /3 = l,m. 

The relations (12.21) or (12.31) are called the complete integrability conditions. 
Theorem 2.2 (Frobenius). Let G C ]R m x IR n be an open subset and 
X a :G^R n = Af„,i(R), X a o/c/ass C 1 , Va = 

a) // the conditions (12. 2p are satisfied, then the PDE system (12. ip zs com- 
plete/?/ integrable. 

b) Let D C ]R m fre an open and convex subset and G = D x M. n . Suppose 
that the following condition is fulfilled: 3R > and inere exist the continuous 
functions ip, ij) : D — > [0, oo) sncn £na£ 



||X a (t,x)|| < p(t)||x||+^(t), Vt G £>, Vx G M n , ||x|| > R, Va = l,m. (2.4) 

(For example, if the PDE (12.11) is linear, then the conditions (12 .4p are sat- 
isfied). 

If the complete integrability conditions (12.21) are satisfied, then: V(io, ^o) G 
DxP, 3x : D -> IR n ; x o/ c/ass C 2 , solution of the PDE system (12.11) and 

X(* ) = Xq. 

Theorem 2.3. Let G C IR m x R" fre an open subset and 

X a :G^R n = M nA (M),X a of class C\ Va = T~m. 

Lei £>i,D 2 C R m be open subsets and y : D 1 ^ W 1 , z : D 2 ^ R n be 
solutions of the PDE system (12. ip . If Did D 2 is connected and there exists 
t G Di n D 2 snca t/iai y(t ) = z(t ), then y(t) = z(t), \/t e D 1 C\ D 2 . 
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Definition 2.2. Let D C R m be an open subset and P a : D — >■ A^ n ,fc 
6e functions of class C 1 . We say that the curvilinear integral / P a (t)dt a 



is path independent (on D), if for any two points t ,ti e D and any two 
piecewise C 1 curves r\ : [a, b] — > D, A : [c, d] — >■ _D ; with 77(a) = A(c) = to 
77(6) = A(d) = ti, we have 



J P a (t)dt a = J P a (t)dt a . 



Theorem 2.4. Let D C M m fre an open subset and P a : D — >■ .M n ,fc(lL) 
C 1 functions, Va = l,m. If D is a convex set (sufficiently, connected and 
simply connected) , then the following statements are equivalent: 

ii) 3^ : D Mnfiffi) solution of the PDE system 

^-(t) = p a (t), Vie A Va = T^. 

raj T/ie curvilinear integral 

' ' P a (t)dt a 



I 



7 

is path independent on the set D. 



In the conditions 1) - Hi), we have: 

a) If £ is a solution of the PDE system of H), and 7 : [a, b] — > D is a 
piecewise C l curve, then 



/ 



P a (t)dr = aj(b))-^(a)). 



b) Let t G D be a fixed point. For t G D, let 7 t0jt : [a, 6] — >■ D be a 
piecewise C 1 curve, from rj to t. The primitive 

A* n , fe (K), £(*)= 1 P a (s)ds a 

7t 0l * 

is a solution of the PDE system of n). 
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For example, if D is a star-shaped set with respect to £ > then the primitive 
£ can be written alternatively 

i 

£(t) = J{r- t%)P a ((l - r)to + rt) dr. 
o 

In case that the PDE system (12.1 j) is linear, i.e., 

— = M a (t)x + F a (t), Va = l,m, (2.5) 

with M a : D ->■ A4„(R) and F a : D ->■ M n = M nA (R) of class C 1 , the 
complete integrability conditions become 

-^(t)x + °-^{t) + M a (t)(M,(t)x + F p (t)) 
= -^{t)x + ^f(t) + Mp{t){M a (t)x + F a {t)), 



Vt e D, Vx E W 1 , Va,/3 = l,m, 

which is equivalent to 



^(t) + M a {t)Mp{t) = ^(t) + M p (t)M a {t) (2.6) 



+ (*) = M p (t)F a {t) + ^f(t). (2.7) 

We have obtained the following result: 

Theorem 2.5. Let D C ]R m fre an open and convex subset, let M a : D — > 
M n (M) be C 1 matrix functions, Va = l,m and let F a : D — > R n = M n ,iQ&) 
beC 1 vector functions, Va = l,m. Suppose the relations (12. 6p . (12. 7p are true, 
V£ G -D, Va, /3 = 1, m. Then the problem 

dx 



— = M a (t)x + F a (t), Va = l,m, 

x(t ) = Xq 

has a unique solution x : D — > M. n . This solution is of class C 2 . 
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Further, everywhere, D will be an open and convex subset of R m , and 
M a : D — > Ai n (M), Vet = l,m, are matrix functions of class C 1 , which 
verifies the relations (12.61) . Wt G D, Wa, {3 — 1, m. 

There exists a unique matrix solution 

X (-,t ):D^M n (R) 

of the problem 

— = M a (t)X, Wa = l,m (2.8) 

X(t ) = I n . 

(For those n problems equivalent to the matrix problem, we apply the The- 
orem I2.5p . 

Definition 2.3. The matrix function 

x (.,.):DxD^M n {R) 
is called the fundamental matrix. 

Proposition 2.1. The fundamental matrix has the following properties: 

a) x{t,ta)x{toM) = x(t,ti), V*o,ti,t £ A 

b) x(to,t ) = I n , Wt G D, 

c) xit,^)" 1 = x{to,t), Wt ,teD. 

d ) ^(x(to,t)) = - X (to,t)M a (t), Wt G D, Wa. 
Proof, a) UY(t) = x(Mo)x(*o, then 

ay d 

= ^(x(t,t ))x(to,ti) = M a (t) X (t, t )x(*o,ti) = M a (t)F; 

y(*o) = x(*o,*o)x(*o»*i) = 4x(*o,*i) = x(*o,*i). 
Hence y(i) and %(t, £i) are both solutions of the matrix PDE system 

— = M a {t)X, Va = l,m, 

which coincide for t = t . From uniqueness it follows that Y(t) = x(t, h), Wt. 

b) Direct consequence of the definition of the function x(£,t )- 

c) It follows readily from a) and b). For a), we take t\ = t, etc. 

d) Differentiating the identity 

X(Mo)x(*0,*) = In 
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with respect to t a , we find 

■^;{x(t,to))x{to,t) + x(t,t ) — ( X (to,t)) = 

or 

M a (t) X (t, t )x(t , t) + X (t, to)-^(x(to, t)) = 0, 

i.e, 

X(t,t )-^(x(t ,t)) = -M a (t). 
Multiplying at the left-hand side by x(Mo) 1 = x(^o^)) we get 



d -(x(t ,t)) = -x(to,t)M a (t). 



dt 



Proposition 2.2. The Cauchy problem 



□ 



— = M a (t)x, Va = l,m, (2.9) 
x(t Q ) = x 

has the solution x : D — > W 1 , x(t) = x{t, to)xo- 

Definition 2.4. Let us consider the PDE system (12.91) . The homogeneous 
PDE system 

^-(t) = -Ml{t)y(t), Va = lT^ (2.10) 
is called the adjoint system. 

The complete integrability conditions of the adjoint system are 

a + MlMJ = -^J- + Ml Ml 



~~ dtP a p ' " " dt a p a 
or 

dM « + M,M a = - d -^ + M a M, 



dtfi p dt< 
dM a %r xe 8Mb 



+ M a M p = — f + M p M a , 



i.e., identical to the relations (12. 6p of complete integrability of the system 

flZSD. 
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Proposition 2.3. a) The matrix solution of the Cauchy problem 

— = -MjX, Va = l,m 
X(t ) = I n 

is <5>{t,t ) = x{to,t) T . 

b) The solution of the adjoint Cauchy problem 

<p(to) = ifo 

is 

<p(t) = $(t,t )<po = X (to,t) T ¥o- 



Proof, a) We use the Proposition 12. 1\ d), i.e., 

-^(x(to,t)) = - X (to,t)M a , 

which is equivalent to 

^(x(to,t) T ) = -M2x(to,t) T . 
b) follows immediately from a). □ 



Theorem 2.6. In the conditions of Theorem \2.5\ the solution of the Cauchy 
problem 

— = M a (t)x + F a (t), Va = l,m, (2.11) 
x(t ) = x 

is 

x:D^R n , x(t) = x(t,t )x + J X (t,s)F a (s)ds a , 

7t ,* 

where jt ,t is a piecewise C 1 curve, included in D, covered from to to t. 
The curvilinear integral / s)F a (s)ds a is path independent. 
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Proof. We show that the curvilinear integral is path independent. According 
to the Theorem 12.4} we must show that 

or 

OF OF 

- X (t, s)Mp(s)F a (s) + X (t, s)q^(s) = ~x(t, s)M a (s)Fp(s) + X (t, s )g^( s )> 

and these are equivalent to the relations (12. 7p . 
Now we write the sheet x(t) as 



x(t) = x(t, t )x + x(t, t ) J X (to,s)F a (s)ds a . 

n ,t 

According to the Theorem 12.41 we get 

— ( J x{to,s)F a {s)As a ) = X (to,t)Fp(t). 



it Q ,t 

It follows that 



Ox 

= Mp(t) X (t,t )X + 



+Mp(t) X (t,t ) J X (to,s)F a (s)ds a + x(t,to)x(t ,t)Fp(t) 

7t .t 



= M p (t) [x(t, t )x + j X (t, s)F a (s)ds a j + F p {t) = M p (t)x(t) + F p (t). 

7t ,t 

One verifies easily the initial condition x(to) = Xq. □ 

3 Controlled nonautonomous linear PDE 
system of first order 

Our main results include generalizations to multitime case of the single-time 
control (see, for example, [2], [9]) in the vision of Lawrence C. Evans, Lev S. 
Pontryagin. They are complementary to the results in [3], [I], [7], [8], [10] 
j. Related topics can be found in the papers PQ, [5]. 
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Let D C M m be an open and convex subset, let M a : D — > Ai n {R) be 
C l quadratic matrix functions, let N a : D — > Ai rit k(M) be C 1 rectangular 
matrix functions, and let u a : D — > M. k = Aik,i(M) be C l vector functions, all 
indexed after a = 1, m. 

We consider the evolution PDE system 



— = M a (t)x + N a (t)u a (t), Va = l,m. (3.1) 
Its complete integrability conditions are equivalent to 



9 -^{t) + M a (t)Mp(t) = d -^(t) + M,(t)M a (t), 

8N du 
M a {t)N,{t)u,{t) + -^(t)u a (t) + N a (t)^(t) 

= M,(t)N a (t)u a (t) + -^(t)u p (t) + Np(t)-£{t), 



Vt G D, Va,/3 = l,m. 

Definition 3.1. Suppose that the matrix functions M a (-) verify the relations 
(12. 6p . Vt G £>, Va,/3 = l,m. TTie vector space 



U 



iu = (w a ) a= i^ u Q : D M fc = A^jfc,i(M), of class C\Va — l,m 

and which verify the relations (13. 2p for all a, /?| 

zs called the control space. 

From the Theorem 12.61 we obtain immediately 

Theorem 3.1. If the matrix functions M a {-) verify the relations (12. 6p . Vt G 
D, Va, /3 = l,m and u = (u a ) a= T^ is a control, then the Cauchy problem 

[)x 

— = M a (t)x + N a (t)u a (t), Va = l,m. 

x(t ) = x (t G D, x G M n ) 

/ias a unique solution 

x:D^R n , x(t) = X (t,t )x + J X (t,s)N a (s)u a (s)ds a , 

7t , i 
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where j to>t is a piecewise C 1 curve, included in D, covered from t to t. 

The curvilinear integral J x{t, s) N a (s)u a (s)) ds a is path independent and 

7 

the solution x(-) is of class C 2 . 

Further, in this paper, D will be an open and convex subset of M. m , the 
C 1 quadratic matrix functions M a : D — > Ai n (M), Va = l,m will verify the 
relations (12. 6p . Wt G D, Va,/3 = l,m and the rectangular matrix functions 
N a : D -> M n>k {R), will be of class C 1 , Va = T^rh. 

Definition 3.2. The pair (s,y), s G D, y G W 1 is called phase of the PDE 
system (13. ip . 

a) Let (to, Xq), (s, y) G DxM". We say that the phase (t , xq) transfers to 
the phase (s,y) if the Cauchy problems { (I3.ip . x(U) = x } and {( 13~Tj) . x(s) = 
y} have the same solution (for the same control u(-)); or, equivalently, the 
solution x(t) of the Cauchy problem { (I3.ip . x(tn) = x } verifies also the 
condition x(s) = y. We will say that the control u(-) transfers the phase 
(t ,x ) into the phase (s,y). 

b) The phase (t,x) is called reachable (respectively pseudo-reachable) if 
there exists a point t G D, with t$ < t a , Va, (respectively, if there exists 
a point t G D, to ^ t), and there exists a control u(-) which transfers the 
phase (to,0) into the phase (t,x). 

c) The phase (t,x) is called controllable (respectively, pseudo-controllable) 
if there exists a point s G D, with s a > t a , Va, (respectively, if there exists a 
point s G D, s ^ t), and a control u(-) which transfers the phase (t,x) into 
the phase (s, 0). 

d) Let to,t G D, with t^ < t a , Wa, (respectively, let to,t G D, t ^ t). 

The PDE system (13. ip is called completely reachable (respectively com- 
pletely pseudo-reachable) from to to t if for any point x G M. n , the phase 
(t ,0) transfers to the phase (t,x), i.e., for any x, the phase (t,x) is reach- 
able (respectively, pseudo-reachable) with the same to- 

e) Let t G D. The PDE system (13. ip is called completely reachable (re- 
spectively, completely pseudo-reachable) at the moment t, if for any t G D, 
with £q < t a , Va, (respectively Vt G D, t ^ t), and for any x G M. n , the 
phase (t ,0) transfers into the phase (t,x). 

f) Let t ,t G D, with t^ < t a , Va (respectively, let t ,t G D, t ^ t). 

The PDE system (13.11) is called completely controllable (respectively, com- 
pletely pseudo-controllable) from to to t if for any point x G M n , the phase 
(t ,x) transfers into the phase (£,0), i.e., for any point x the phase (t ,x) is 
controllable (respectively, pseudo-controllable) with the same t. 
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g) Let to G D. The PDE system ( 13.1 \ is called completely controllable 
(respectively, completely pseudo-controllable) at the moment t , if\/t G D, 
with t a > Va, (respectively, Wt G D, t 7^ t), and for any point x G W 1 , 
the phase (to,x) transfers into the phase (t,0). 

h) The PDE system (13. ip is called completely reachable (respectively, com- 
pletely pseudo-reachable) if it is completely reachable (respectively, completely 
pseudo-reachable) at any moment of D. 

The PDE system (13.11) is called completely controllable (respectively, com- 
pletely pseudo-controllable) if it is completely controllable (respectively, com- 
pletely pseudo-controllable) at any moment of D. 

The multitime control property does not only depend on the dimensions 
m and n but on how matrices M a and N a interact. 

The phase (t ,x ) transfers into the phase (t\,y) -<=>- 3u(-) = (u a (-)) a 
control such that the solution x(-) of the problem { (I3.ip . x(t ) = x } verifies 
also x(ti) = y, equivalent to 

3u(-) = (tt a (-)) a control such that 



x(t) = x(t,t )x + J x(t,to)x(to, s)N a (s)u a (s) ds a and x(h) = y 

Ti ,t 

<^=^> 3m( - ) = (u a (-)) a control such that 
V = x{ti,to)(xo+ J x(to,s)N a ( y s)u a (s)ds a 

T*o.«i 

<^=^> 3m(-) = (« Q (-)) a control such that 
X{t , h)y - Xo = J x{h-,s)N a (s)u a {s)ds a . 

We introduce the set 



V(t Q ,t):={ j X {to,s)N a (s)u a {s)ds c 



(w a )Q, = T^: is a control 



The set V(t ,t) is a vector subspace of M n . It is called the controllability 
space. Since the curvilinear integral is path independent, we remark that 
V(t ,t) does not depend on the curve jt ,t, which joins t to t, but depends 
on the multitimes to and t. Also t )V(t , t) = V(t,t ). 
From the foregoing arguments, it follows immediately 
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Theorem 3.2. Let us consider the system (13.1 1) , with the matrix functions 
M a (-) verifying the relations (12.61) . 

i) The control {u a ) a= i^ transfers the phase (to, ^o) to the phase (t,y) if 
and only if 

x (t ,t)y-x = J x(t ,s)N a (s)u a (s)ds a . 

7t , * 

ii) The control (u a ) a= i^ transfers the phase (t ,x ) to the phase (t, 0) if 
and only if 

Xo+ J x(to,s)N a (s)u a (s)ds a = 0. 

7t , * 

Hi) The phase (to>£o) transfers into the phase (t,y) if and only if 

xq - x(ta,t)y G V(t ,t) 

equivalent to 

y - x(t,t )x G V(t,t ). 

iv) The phase (to^o) controllable (respectively, pseudo-controllable) if 
and only if3t G D, with t a > t$, Va (respectively, 3t G D, t ^t ) such that 

x G V(t Q ,t). 

v) The phase {t,y) is reachable (respectively, pseudo-reachable) if and only 
if 3t G D, with < t a , Vet (respectively, 3t G D, t ^ t) such that 

yeV{t,t ). 

vi) Let to,t G D, with £q < t a , Va (respectively, let t ,t G D, t ^ t). 
The PDE system is completely controllable (respectively, completely pseudo- 
controllable) from the multitime t into the multitime t if and only if 

V(t ,t) = R n , equality equivalent to V(t,t ) = R n . 

vii) Let t ,t G D, with tg < t a , Wa (respectively, let t ,t G D, t ^ t). 
The PDE system is completely reachable (respectively, completely pseudo- 
reachable) from the multitime t into the multitime t if and only if 

V(Mo) = R n , equality equivalent to V(t ,t) = R n . 

viii) Let t ,t G D, with t$ < t a , Wa (respectively, let t ,t G D, t ^ t). 
The PDE system is completely controllable (respectively, completely pseudo- 
controllable) from the multitime to wiio the multitime t if and only if it is 
completely reachable (respectively, completely pseudo-reachable) from t tot. 
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According to the Theorem I2.4[ the curvilinear integral 

x(t ,s)N a (s)Nj(s) X (t ,s) T ds a 
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is path independent if and only if, for any a, (3 = l,m, the following condi- 
tions are satisfied: 

d dN 

( X (t , s))N a Nj X {to, s) T + X (to, s)—^-N2x{to, s) T + 



or 



8N T d 
+X(t , s)N a -^ X {t , s) T + X (to, s)N a Nj — {x{t , s) T ) 

d dN 
= Q^(x(to, s^NpNjxito, s) T + X (to, s)^Nj X (to, s) T + 

+X(t , s)Np-^ X (-to, s) T + X (to, s)NpNj—(x(to, s) T ) 

- X (t , s)MpN a Nj X {t , s) T + x(*o, s)-^-Nj X {to, s) T + 
8N T 

+ X (to, s)N a —f X {to, s) T - X (to, s)N a NjMj X {to, s) T 



dsP 

- X {to, s)M a NpNj X (to, s) T + X (t , s)-^Nj X (to, s) T + 



dNJ T T T 

+X(t , s)Np-^ X (to, s) T - x(to, s)NpN]Mlx{h, s) T . 

Since the fundamental matrix x{to,s) is invertible, the foregoing equality is 
equivalent to 

—MpNaNj + —^Nj + N a —^- - N a NjMj 

OSp OS/3 

T 8N B T dNJ T T 

= -M a N p Nj + -^Nj + N p -^- - NpNjMj. 

In this way, we have proved 

Proposition 3.1. Let to G D, fixed. The curvilinear integral 

X (to,s)N a (s)N ( l(s)x(to,s) T ds a 

7 
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is path independent on D (in the sense of definition 12.2)) if and only if, for 
any a, (3 = l,m, the relations 



dN dN T 
M a NpNj + -^N T a + N a -^- + NpNjMj 

r dN B T dNj T T 



(3.3) 



are verified on D. This is equivalent to, 

dN n „ -r / „ „ T cW, 



(3.4) 



or 



aN / r)N \ T 

T 9N ( T 9N J\ T 



(3.5) 



It is sufficient, for example, that for any a, ft = 1, m, to have 

M*N P NJ + 8 -^Nl = MpN a Nl + ^Nj (3.6) 



or 



M a NpNj + = M p N a Nl + Np^. (3-7) 

Proposition 3.2. Let us suppose that the matrices M a (-) verify the relations 
flXHj) . \/t E D, Va, P = T~m. We fix t E D. For each v G R™ and a = T~m, 

we consider the functions 

u a>v : D ->• Ua>w (s) = N^(s) X (t , s) T v, Vs G D. 

The following statements are equivalent 

i) For any v E MJ 1 , the family (u a)V ) a= Y^ is a control for the PDE system 

m>- 

ii) For any a,/3 = l,m, the relations (13. 3 p are satisfied on the set D, i.e., 



M a NpNj + -^Nl + N a -^ + NpNjMl = 

T dN g T ON] 
= M p N a Nl + -^Nj + Np-^. + Na NjMj. 
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Hi) The curvilinear integral 

x(t 0l s)N a (s)N2(s) X (to,s) T ds a 



is path independent on the set D. 

Proof. The family (u a ,v) a =n^, is a control if and only if it verifies, Va, j3 
l,m, the relations (13. 2p on the set D, i.e., 



M a N p N p x (to, s) l v+ ^-Nj X (to, s) T v+ 



M p N a Nlx{h, s) T v + -g^Njxita, s) T v+ 



+N a ^ X {to, s) T v + N a N> (x(*o, s) 1 )v 

.a P 

+ N ^x(to, s) T v + N Nj-^( X (t o , s) T )v 

equivalent to 

M a NpNjx{to, s) T v + -^-Nj x (to, s) T v+ 

+N a ^x(to, s) T v - N a N2Mj X (to, s) T v = 
= M p N a NZx{tv, s) T v + —JLNjxfa, s) T v+ 



ds c 

ONJ 
ds c 



+ N P^x(t , s) T v - NpNjMT x (t , s ) T v 



or 



dN dN T \ 

M a N p Nj + -^Nl + N Q ^ - N a N^Mj) X (to, s) T v = 

f t 3N R -r ONJ -r t\ 

[M,N a Nj + -^Nj + N p -^ - NpNjM^xito, s) T v, 
equivalent to 



/ dN dN T \ 

\M a NpNp + -^jNl + N a -^ + NpNjMjj X (to, s) T v 

{ t 8Nr -r 9NJ -r -r\ 

[MpN a Nl + -^Nj + N p -^ + N a NjMj) X (to,s) T v. 



(*) 
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The implication ii) ==>- i) follows immediately. 

Let us prove i) =>■ ii). Since for each v G W 1 , the family {u a>v ) a= Y^i is 
a control, it follows that the relations (*) hold for any v G M n , whence we 
deduce that for any matrix A G JA ntP (R), Vp, we have 

\M a N p Nj + -^-Nj + N a -^ + N p NjMZ)x(to> s) T A 
Taking A = I n , we find 



-r dN B ^ ONl -p -p\ 



/ 8N 8N T \ 

(M a NpNj + -^N T a + N a -^- + NpNjMtyfo, s) T 

f t ONr t dNl _ T \ 



The matrix x(to> S ) T is invertible. The last equality is multiplied in the 
right-hand side by (x(to, S ) T ) _1 ; obtaining the relation (13. 3p . 

The equivalence of the statements ii) and Hi) is just the Proposition 

i3~n □ 



Definition 3.3. Suppose that, for any a,/3 — l,m, the relations (I3.3P are 
true. The matrix function 



C:DxD^M n (R), C(t ,t):= J X (h, s)N a (s)N>(s) x (t , s) [ ds a 

is called the controllability gramian. 
The matrix function 

K(t ,t):= J x(t,s)N a (s)N2(s) X (t,s) T ds a 

7*o,< 

is called the reachability gramian. 

The controllability gramian is used to determine whether or not a linear 
PDE system is controllable. The reachability gramian is used to determine 
whether or not a linear PDE system is reachable. One observes immediately 
that 

TZ(t ,t) = - C(t,t ), Vt ,teD 

and 

X (t, t )C(t , t) X (t, t ) T = - C(t, t ), Vt , t G D. 
Hence the matrices C(to,t), C(t,to), lZ(to,i), TZ(t,t ) have all the same rank. 
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Definition 3.4. Let A G A4 Ptq (M.) be a real matrix. Denote lm(A) and 
Ker(^4) ; the image, respectively the kernel of the linear map 

f : R q = M qi i(R) ^W p = M p> i(R), f(x) = Ax. 

Of course, the subset Im(A) is a vector subspace of Ai Pt i(M) generated by the 
columns of the matrix A. 

Theorem 3.3. In the conditions of Theorem \'S.1\ if\/a,/3 = l,m, the con- 
ditions ( 13. 3 p are true, then, for any t and t with t a > t%, Va = l,m (or 
t a < £q ; Va = l,m), we have 

V(t ,t) = lm(C(t ,t)). 

Proof. The inclusion V(to,t) C lm(C(to,t)) is equivalent to 

(V(t ,t)) ± 2 (Im(C(t ,t))) ± = Ker((C(to,t)) T ). 

We have b G Ker ((C(t ,t)) T ) <^=> ((C(t ,t)) T )b = <J=> 6 T C(t ,t) = 0. 
Hence b T C(t ,t)b = 0. 

The controllability gramian is independent on the curve 7 covered from 
the multitime to to the multitime t. Particularly, we fix 7 as being the straight 
line segment which joins the points to, t, i.e., 7(1") = r(t — 1 ) + t , r G [0, 1]. 
It follows 

\ m 

C(t ,t)= / ^(t Q -^) X (to,7(r))iV a (7(r))iV Q T (7(r))x(to,7(r)) T dr 
a=1 

1 

= / E r-^o)x(to,7(r))iV a (7(r))iV Q T (7(r))x(to,7(r)) T dr 

a With *V*o 

On the other hand we get 
1 

J E ^ ~ ^o)^x(to,7(r))iV Q (7(r))iV a T (7(T))x(to,7(r)) T &dr = 0, 
a with t Q >^ 

or 

1 

/ E (* a - 

a with t Q >tg ! 
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^0,7(^)^(7^)) 



dr = 0. 



It follows that for any a with t a ^ t% and, Vr e [0, 1], we have 



6 T X (t ,7(r))iV Q (7(r)) = 0. 



Hence 



bT / ^x(to,7(0)^(7(r))n Q ( 7 (r))7 a (r)dr 



1 




o a with t a ^t; 



Consequently b G (V(to ) ^)) ± - 

Let us prove the inclusion lm(C(t Q ,t)) C V(i ,t). For that, we select 



w = C(t 0j t)ueIm(C(i ,*)), u> = / x(io,s)^ Q (s)^(s)x(to,s)%ds Q . 



According to the Proposition 13. 2\ the family (« a ) a= r^ is a control, hence 



If, for any a,/3 = l,m, the relations (13. 3p are true (hence we have also 
V(t ,t) = lm(C(t ,t)) for t a > (or Va), then from the Theorem O it 

follows 

Theorem 3.4. Suppose that we are in the conditions of the Theorem \'d.l\ 
and furthermore, for any a,/3 — l,m, the conditions (13. 3 j) are true. 

i) Le£ t ,t e D, t ^ t, 1% < t a , Va or t a < *g>, Va. T/jen toe p/jase 
(to,^o) transfers to the phase (t,y) if and only if 




We choose 



u a (s) = N a (s)x(t ,s) v. 



w = 




□ 



- x{to,t)y e lm(C(t ,t)) 



which is equivalent to 



y - x{t,to)xo e lm(C(t,t )). 
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ii) The phase (t Q , x ) is controllable if and only if 3t G D, with t a > tifi, 
Wa, such that 

x G lm(C(t ,t)). 

Hi) The phase (t,y) is reachable if and only ifBto G D, with t$ < t a , Va ; 
such that 

yelm(C(t,t )). 

iv) Let t ,t G D, with t% < t a , Va. The PDE system is completely 
controllable from to to t if and only if 

rank C(t , t) = n ( W = Im(C(t , t)) 

v) Let t Q ,t G D, with £g < t a , Va. The PDE system is completely reach- 
able from t into t if and only if 

rank C{t ,t) = n ( R n = Im(C(* , *)) )■ 

4 Controlled autonomous linear PDE 
system of first order 

A very special case is those of controlled autonomous linear PDE system 
of first order, when the matrix functions M Q , N a are constants. Then the 
relation (12.61) becomes 

M a M p = M ;.!/„. 

the relation (13.21) reduces to 

M a N p u p + N Q ^ = M p N a u a + 

and the relation (13. 3p can be written as 

M a NpNj + NpNjMj = MpN a Nl + N a NjMj. 

On the other hand, we have 

x (t,t ) = e M ^ ta - t o\ 

In this case, the fundamental matrix function to) is defined for any 
(Mo) G R m x R m , and not only for (t,t ) G D x D. Moreover, x(Mo) = 
X (t-t ,0). 
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The controllability gramian matrix becomes 

C(to,t)= [ e- M ^-^N a Nl(e- M ^- tP A T ds^ 



7i ,t 



Since the relations (I3.3P are verified on M. m , from the Proposition 13 . 1 1 it follows 
that the curvilinear integral is path independent on R m . Hence the control- 
lability gramian C(t , t) is defined on R m x R m and also C(t , t) = C(0,t-t ). 
In fact, if the matrix functions M a and N a are constant Va, one can take 



D = R m . 



Definition 4.1. For a = l,m, let us consider the constant matrices M a e 
M n iR), N a G M n ,kQ&), such that 



M a Mg = MpM a , Va, /? = 1, m. 



For each a — 1, m, we define the matrix 

G a = [N a M 1 N a M 2 N a . . . M m N a . . . Mf 1 ■ ■ . . . ■ M%> ■ N a ... 

made from all block matrices of the form 

Mf 1 ■ M* 2 ■ ■ N a 

with < hi, k 2 ; ■ ■ ■ ; k m < n — 1. Further, we need to specify the order in 
which one ranges the block matrices M-f 1 ■ M 2 2 ■ . . . ■ M^ m ■ N a in the matrix 
G a . In this way, the matrix G a will be well defined. 
On the set 

{Ox; k 2] . . . ; km) e N m | < h; k 2] . . . ; k m < n - l}, 
we define the following order relation, denoted by ^: 

(h; k 2 ; ■ ■ ■ ; k m ) di (qi; Q2', ■ ■ ■ ; q m ) if 



h + k 2 ^ Vk m <q 1 + q 2 J ^ h q m 

or 

h + k 2 -\ h k m = qi + q 2 H h q m and 



k\ > qi or 3j — 2, m such that k\ = q±, k 2 = q 2 , . . . , kj-i = qj-i, kj > qj 



or 
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{ki,k 2 ; ...;k m ) = (q 1 ;q 2 ; ■ ■■;q m ). 
One verifies quickly that ^ is an order relation. The block matrices 

M* 1 ■ M 2 fc2 • ■ ■ ■ N a 

are ranged in G a in increasing order of (k\\ k 2 ; . . . ; k m ), relative to the order 
relation ^. This means in fact that the block matrices are written in the 
increasing order of the sum ki + k 2 + • • • + k m ; in case that two such sums 
are equal, the block matrices are written in lexicographic decreasing order of 
\k\] k 2 ; . . . ; k m ) . 
The matrix 

G \— ^ G\ G 2 ... G m J . 
is called the controllability matrix of the PDE system (13. ip . 
Theorem 4.1. Let t ,t G W 71 such that t% <t a ,Va = T~m. Then 

rank C(t, to) = rank C(to,t) = rank G, 
relation equivalent to 

Im C(t ,t) = Im G, 

or 

Ker(C(t ,t) T ) = Ker(G T ). 

Proof. We already have seen that rank C(t,t ) = rank C(ta,t). Hence, it is 
enough to show the equality Ker(C(t , t) T ) = Ker(G T ). 

The inclusion Ker(C(t 0l t) T ) C Ker(G T ): let b such that (C(t , t)) T b = 0, 
or b T C(to,t) = 0. Consequently we have and b T C(t ,t)b = 0. 

Let s such as £g < s a < t a , Va = l,m, s fixed but arbitrarily chosen. 
Instead of the curve j to ,t we select the union of the segments [to, s] and [s, t], 
covered from to to s, respectively from s to t. Let us parameterize these 
segments by 

Ai(r) = (l-r)to + rs, Vr G [0,1], 
X 2 (t) = (l-r)s + Tt, VrG [0,1]. 
Obviously 7 t0jt is a piecewise C 1 curve. From 6 T C(t , t)b = 0, we find 
i 

J b T e- M ^) ~ OjVJVj ( e -^(A?(r) - ^)) T fe (a « _ ^ dr 
o 

+ / b^e- M ^{r) - O^^e"^ (r) - $) T b {t « _ s «) dr = 0; 
o 
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b T e -M^(T)-4) 



+ 



b T e -MpM{T)-1$) Na 

Since the integrand is positive, it follows 



(t a - s a ] 



dr = 0. 



a=l 



3 f \ 8\ 2 

+ ^ 6 T e- M ^ A 2( T ) -^o)^ (t a -s a ) = 0, Vre[0;l]. 



But s a - 1% > and t Q - s Q > 0, hence 



fe T e -ik^(A 1 (r) - t ) Na = 0> Vr e [ . ^ V a = 1, 



and 



6T e -M /9 (A 2 (r)-i ) JVa = o ) VrG[0;l], Va = l,m. 
We set r = 1, hence Af(l) = s 13 and one obtains 



6'e" 



l,m; 



valid equalities W s 13 G (t^;^), /3 = l,m. But, taking into account the 
continuity, it follows 



6'e" 



-M p {sP-i) N 



a = 0, Va = 1, m, 



(4.1) 



and Vs such that 4 < s 13 < , V/3. 

Let < ki, hi-, ■ ■ ■ , k m < n — 1. In (14.1 1) we differentiate with respect to s 1 of 
k\ times, with respect to s 2 of k 2 times, . . . , with respect to s m of k m times 
(we take also into account the relation (12.61) ): 

(_ 1 )fci+fc 2 +...+fc m6 T e -M /3 (/ - 4) M k 1M k 2 M^N a = 0, 
while for s = to we have 

b T M^M^ ■ ■ ■ M^ m N a = 0, 
hence also b T G = or G T b = 0, i.e., b e Ker(G T ). 
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The inclusion Ker(G T ) C Ker (C(t ,t) T ): let b such that G T b = 0, or 
b T G = 0. It follows that 

b T M kl M k2 ■ ■ ■ M^N a = 0, V0 < k u fa, . . . , k m < n - 1. 

From the Hamilton-Cayley Theorem, and taking into account the relations 
M a M{s \ 1 1 \ /„ . we deduce 

fM^Mj? ■ ■ ■ MtN a = 0, V A*, fa, . . . , k m > 0, 

6 T (^ - s^M^itl - s 2 ) k2 M% 2 • • • (C - s m ) km M^N a = 0, 



valid equalities Vfci, &2, • • • , k m > 0. Hence we have Vp > 0, Va = 1, m, 

m 

6 T (^Af^(^-/)) P iV Q = 0. 



/3=1 

fo T e M^(*5 - 8P) Na = J- ±b T (Mp(4 - s^) P N Q = 0. 

p=0 ^ ' 

It follows that b T C{t , t) = or (C(t , t)) T b = 0, hence 6 e Ker(C(t , t) T )- □ 
Remark 4.1. From the proof of the Theorem \4.1l we see that the inclusion 

Ker(G T ) C Ker(C(t , t) T ), i.e., Im C(t ,t) C Im G, 

is true /or any £ hence we have also 

rank C (t , < rankG, V* 0) * € M m 

Example 4.1. There exist linear autonomous PDE systems for which we 
have to,t G M. m , to ^ t, t% < t a , Va, such that rankC(t ;£) < rank G. 
Indeed, let us take 

m = 2, n = 2, k = l, D = R 2 

These verify M l M 2 = M 2 M 1 = and AfiiVi = JVi, M l N 2 = 0. 
Control space: u = (u±, u 2 ) is a control if and only if 

M x N 2 u 2 + iVx-^ = M 2 N lUl + N 2 ^; 

Ot 2 Oti 
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the relations MiN 2 = 0, M 2 = transforms the foregoing PDE in 

ot 2 Oti 

equivalent to 

du 1= du 2 = 

dt 2 ' dt 1 

Consequently u = (ui,u 2 ) is a control if and only if there exist f\, f 2 
R — >> R, of class C 1 such that 

uiCt 1 ,* 2 ) = /iCt 1 ), « 2 (tV 2 ) = / 2 (t 2 ), vet 1 ,* 8 ) e m 2 . 

The relation 

MiN 2 N 2 + N 2 Nj Mj = M 2 N x Nj + N x Nj M 2 T 

(see the condition (13.31) ) is obvious since M 1 N 2 = and M 2 = 0. 
The rank of the matrix 

G ( .V; MiiVi .U>.Y, MiM 2 Ni N 2 M±N 2 M 2 N 2 M X M 2 N 2 ) 

1 1 
v 1 

is 2. 

We compute the matrix C(t , *), with t = = (0, 0), t = (t 1 , 0), (t 1 ^ 0), 
i.e., tg = t 2 = 0. For that, we select 7(1") = (r, 0), r G [O,* 1 ]; 7 is a curve 
joining the two-time (0, 0) with t = (t 1 , 0). Then 

i 1 

C((0,0);(t\0)) = / e" rMl • ^N? ■ (e~ TMl Y dr. 





But 

C((0, 0); (*\ 0)) = I ( 6 7 ) ( e~ T ) dr = | ( ^ ° ) dr; 



C((0,0);(t\0)) 



16 

2 u 





Consequently the rank of the matrix C((0, 0); (t 1 , 0)) is 1, strictly smaller 
than the rank of G. 

The Theorem 13.41 can be rewritten as 
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Theorem 4.2. For a = l,m, let us consider the constant matrices M a G 
M n (l), iV a G M n<k (R), such that 



M a Mfi = MpMu, Va, p = l,m, 

M a NpNj + NpNjMj = M?N a Nl + N a N^Mj, Va, /? = T^m. 
VFe consider the autonomous PDE system 

dx 



M a x + N a u a (t), Va = 1, ra, 

where u = (u a ) a= ji^ is a control, i.e., u a : W" 1 — > M, k = A^fc,i(M) o/ c/ass 
C 1 , Va, and 

M a ^(t)+^ a ^(t) = M^iV a « Q (t) + ^^(t), VtGK m , Va,/3 = I^n". 

Let G be the controllability matrix of this PDE system. 

i) If t a > (or Va, then the phase (t ,xo) transfers to the phase 
(t, y) if and only if 

Xo _ e M a (t% -t a ) y e Im ( G! ) > 

equivalent to 

y- e M ^ ta -^)x elm(G). 

ii) The phase (to,^o) ^ s controllable if and only if x G Im(G). 

One observes that if exists a multitime to such that the phase (t ,xo) is con- 
trollable, then for any multitime t, the phases (t, xq) are controllable. 

Hi) The phase (t,y) is reachable if and only if y G Im(G). 
One observes that if there exists a multitime t such that the phase (t, y) is 
reachable, then for any multitime s, the phases (s, y) are reachable. 

iv) If the phase (to,Xo) is controllable (or reachable), then for any multi- 
time t, the phases (t, xq) are controllable and reachable. 

v) Let t ,t G D, with < t a , Va. The PDE system is completely con- 
trollable from t to t if and only if rank G = n. 

vi) Let t ,t G D, with < t a , Va. The PDE system is completely 
reachable from t to t if and only if rank G = n. 

vii) If there exist t ,t G D, with < t a , Va and if the PDE system 
is completely controllable (or completely reachable) from to to t, then the 
PDE system is completely controllable and completely reachable (equivalent 
to rank G = n). 
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Example 4.2. Let us give an example of multitime linear PDE system, with 
constant matrices M a , N a , for which we have rank G = n, but no state 
(to,xo), with xq 7^ 0, is controllable. Let us consider 



Mi = M 2 = M 3 n = M 






















V o l 







D 


= I 


i 3 ; 








l 


1 











1 







N 3 






The foregoing matrices verify the relations 

MN X = N 2 , MN 2 = N 3 , MN 3 = N x 

and the relations (12.61) are obvious. 

Let us determine the control space: u = (ui,u 2 ,u 3 ) is a control if and 
only if the relations (13.21) are true. If in (13.21) we put a — 2, /3 — 3, then 

M 2 N 3 u 3 {t) + N 2 ^{t) = M 3 N 2 u 2 {t) + N 3 -^{t), Vt G M 3 , 

or 

Mt)Ni + ^jjtQW = Mt)N 3 + ^(t)N 3 , Vt G R 3 . 
Since N x , N 2 , N 3 are linearly independent, it follows that 

u 3 (t) = 0, u 2 {t) + -^(t) = 0, Vt g M 3 , 

hence u 3 (t) = u 2 {t) = 0, Vt G M 3 . 

If in (13.21) we take a = 1, (3 = 2, then we get 

M x N 2 u 2 (t) + N x ^{t) = M 2 N lUl {t) + N 2 -^(t), Vt G M 3 , 

or 

Ni-^(t) = N 2Ul (t), Vt G K 3 , 

hence «i(t) = 0, Vt G IR 3 . We have proved that the control space is {0}, and 
consequently 

V(t ,t) = {0}, V(t ,t) G M 3 x M 3 . 

These results and the Theorem 13.21 show that no phase (t , x ), with x ^ 0, 
is controllable or reachable. 

We remark that the controllability matrix G has the rank 3 = n. The 
conclusion of the Theorem 14.21 is no longer valid. The reason is that the 
conditions ( 13. 3 p are not true. 
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5 Comments 



In this paper, the functions which define the PDE systems (for example ( 12. ip . 
(13. ip . etc.) are of class C 1 and satisfy the complete integrability conditions 
(of type (TjL2D, dZSD + dZZD, (ESD + tD)- Also, the general relations flZH 
are verified - we have linear PDE systems. Hence, throughout, the Cauchy 
problem { (13.11) . x(t ) = x } has a unique solution, global defined, and it is 
of class C 2 (Theorems I2TT1 1^21 123]) . 

Sometimes, in the papers [12] - [23] , the functions which define the PDE 
systems (for example, the controls) are piecewise C l functions; in this case, 
the complete integrability conditions are piecewise satisfied. Identically, the 
solutions will verify the PDEs in the piecewise sense. Generally, the solutions 
are not continuous functions. 

In the paper [13] it is indicated a construction of the solution of a Cauchy 
problem associated to a linear PDE system. It is similar (in a certain sense) 
with those obtained in the Theorem 13.11 But, in this context, the Cauchy 
problem has not a unique solution. To maintain this idea, we give the fol- 
lowing example: let m = 2, n — 1 (hence x(-) = xi(-)), k = 1, D = M 2 , 



M 1 (t) = M 2 {t) = 0, Nx(t) — N 2 (t) — 1, VtG 



o2. 



M*,?) = \l'*t + i-] ; ^ 2 (t\t 2 ) = o,v(t\t 2 )GR 2 . 

Set to — (0, 0), Xq = 0, i.e., x(0, 0) = and we formulate the Cauchy problem 

dx dx 

w = u l (t\t 2 ), — = 0, x(0,0) = 0. 

This PDE system satisfies the piecewise complete integrability conditions 
( 12. 6p and ( 13. 2 p (this can be easily checked); the conditions (13. 2p are true 
there where the control "Ui(-) is of class C 1 , i.e., on the non-connected set 



e\{(t l ,t 2 )lt l +* 2 = i}. 

arbitrary fixed two-time, then 



Here, we have x(t,s) = 1, V(t, s) e M 2 x M 2 . If t = (t l ,t 2 ) e it- is an 



x{t\t 2 ) = X (t,t )x + J X (t,s)N a (s)u a (s)ds a = J Ul (s\s 2 )ds l . 



7*Q, * 7*o, t 



For t 1 + 1 2 < 1, we obtain obviously xit 1 ^ 2 ) = 0. It remains to study the 
case t 1 + 1 2 > 1. 
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Let a G K. be an arbitrary point. We consider the curve 7t 0ji consisting 
in two segments: the first is the segment which joins the point t = (0, 0) 
to the point (a, 1 — a), on the straight line s 1 + s 2 = 1, where wi(-) is 
discontinuous; the second segment joins the point (a, 1 — a) to the point 
(t l ,t 2 ). The first segment, without the point (a, 1 — a), is situated in the 
semiplane {(s 1 , s 2 )^ 1 + s 2 < 1}, and the second, without the point (a, 1 — a), 
is included in the semiplane {(s 1 , s 2 )^ 1 + s 2 > 1}. A parametrization of the 
second segment is 

s^r) = (1 - r)a + rt 1 , s 2 (r) = (1 - r)(l - a) + rt 2 , r E [0, 1]. 

Taking into account that iti(-) vanishes on the first segment, we find 



^(t 1 ,* 2 ) = J \.{t l - a)dr = t 1 
o 

Consequently, the solution is given by the formula 

x{t\t 2 ) -- 

on the set 



o. 



t 1 - a, if t 1 + t 2 > 1 
, if t 1 +t 2 < 1, 



R 2 \{(t 1 ,t 2 )|t 1 + t 2 = 1} 

and x(0, 0) = 0. Here we recognize an infinity of solutions since a is an 
arbitrary point. 

The foregoing solution x(-) can be extended to a continuous function at 
(a, 1 — a), but in rest the function x(-) is discontinuous on the straight line 
{(t 1 ,^)^ 1 + t 2 = 1} (for any given values on this straight line). 
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